Experimental demonstration of a technique to generate arbitrary quantum 

superposition states 
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Using a single, harmonically trapped ^Be^ ion, we experimentally demonstrate a technique for 
generation of arbitrary states of a two-level particle confined by a harmonic potential. Rather than 
engineering a single Hamiltonian that evolves the system to a desired final sate, we implement a 
technique that applies a sequence of simple operations to synthesize the state. 

PACS numbers: 03.67.Lx, 32.80.Qk 



The goal of deterministically synthesizing or "engineer- 
ing" arbitrary states of a quantum system is at the heart 
of such diverse fields as quantum computation [Q and 
reaction control in chemistry For harmonic oscilla- 
tor states, particular non-linear interactions can be used 
to generate special states such as squeezed states. How- 
ever, it is intractable to realize a single interaction re- 
quired to create an arbitrary state. Law and Eberly |^ 
have devised a technique for arbitrary harmonic oscillator 
state generation that couples the oscillator to a two-level 
atomic or "spin" system and applies a sequence of oper- 
ations that use simple interactions. We demonstrate this 
technique on the harmonic motion of a single trapped 
^Be+ ion and include the generation of arbitrary spin- 
oscillator states [Q. Such quantum state control is rele- 
vant to the scheme for constructing a quantum computer 
using trapped atomic ions ||, ||, where we must control 
the quantized micro- mechanical system composed of the 
collective ion normal modes that are used as a data bus 
to transfer information between the ion qubits. These 
techniques could also be used to create input states for 
quantum computing schemes that use continuous vari- 
ables 10] including the code-words that are required for 
fault tolerant computation |^. 

Arbitrary quantum state synthesis is difficult unless 
certain conditions are met. As an example, consider 
a simple quantum system with four energy eigenstates 
labelled |0), |1), |2), and |3). If the system is ini- 
tially prepared in |0), and if couplings that create su- 
perpositions ai\0) + (3i\i) [i — 1,2,3) can selectively be 
turned on, then we can create arbitrary superpositions of 
the form cq |0) + ci |1) -I- C2 |2) -|- C3 |3). Here, the c; are 
complex and subject to the usual normalization condi- 
tion \cif' = 1. This method could be realized in an 
atomic system if the four states were non-degenerate lev- 
els with different energy separations and coherent tran- 
sitions |0) \i) could be driven by applied radiation. 
These requirements are not often met in practice. For 
example, it may be impossible to realize all of the desired 
couphngs |0) Also, if the eigenstates are equally 

spaced like the first four energy levels of a harmonic os- 



cillator, then driving the |0) ^ 
successive transitions |1) ^ |2) 
fixed relations between the &. 



1) transition also induces 
|2) ^ |3), etc. leading to 



It has long been recognized that certain interactions 
can cause harmonic oscillators to evolve to particular 
desired states For example, if the oscillator is ex- 
cited from its ground state with the nonlinear force 
Fqz cos{2Lot) then a "vacuum-squeezed" state is created. 
Such states can be used to increase measurement preci- 
sion in specific applications such as interferometry [ p^ . 
However, it is usually intractable to find the desired force 
or interaction that will create a state with arbitrary co- 
efficients. To circumvent this problem, schemes have 
been proposed [|ll|, |l^ that sequentially couple atomic 
superposition states to the field of a cavity mode and 
statistically prepare arbitrary field states through pro- 
jective measurements. An alternative method has been 
proposed to deterministically map a previously prepared 
superposition of atomic Zeeman states onto the field of 
a cavity [p^ . A more general deterministic scheme to 
prepare arbitrary field states has been suggested by Law 
and Eberly ^, The idea relies on coupling the har- 
monic oscillator to an auxiliary two-level quantum sys- 
tem through a sequence of simple interactions. 

Consider an auxiliary system consisting of two inter- 
nal states of an atom which we label ||) and |t) in anal- 
ogy with the two-level system resulting from a spin-1/2 
magnetic moment in a magnetic field. In practice, the 
harmonic oscillator could correspond to a single mode 
of the radiation field [0 or the mechanical oscillation 
of a trapped atom jisf . The combined energy lev- 
els for this system are depicted in Fig. 1. As an ex- 
ample, we summarize the procedure to create the state 
li) X]i=o...3 I*) starting from the ground state |i) |0). 
It is simplest to first think about solving the inverse 
problem creating the state ||) |0) from the initial 
state II) X]j=o 3 !*)■ procedure begins by apply- 
ing a resonant pulse of radiation that carries out a "tt- 
pulse" II) |3) It) |2) leaving no amplitude in the ||) |3) 
state ( "clearing it out" ) and placing amplitude C|3 in the 
It) 1 2) state. Applying this radiation also causes transi- 
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tions between other states ||) |n) and |t) |n — 1). Because 
in general the coherent transition rates (Rabi rates) are 
not the same for different values of n in this first 
step the amplitudes of the other states change accord- 
ing to Ci„ li) \n) ^ c'^„ li) |n) + d'^^^_, |T) |n - f). The 
Law/Eberly method succeeds because the state | j) |0) re- 
mains unaffected since the state \1) | — 1) does not exist. 

The second step is to induce the transition 
Ci3 IT) |2) + c[., \i) |2) ^ c'/2 li) |2), thereby clearing out 
the It) |2) state. The duration and phase of the sec- 
ond pulse are chosen according to the known values of 
C|3 and c'|2 in order to collapse the superposition state. 
The first two steps have cleared out the ||) |3) and |t) |2) 
states but, in general, non-zero amplitudes remain in the 
|T)|0), |T)|1), |i)|0), U)|f), and |i)|2) states. How- 
ever, by repeating this two-step clearing-out process for 
successively lower values of n, the state amplitudes are 
transferred down the dual ladder of states eventually to 
the ground state |J,) |0). Finally, to achieve the origi- 
nal goal, we apply these same steps in a time-reversed 
fashion to carry out the mapping to |i)X]i=o a'^iiN)- 
In the experiments described below, we demonstrate 
the Law/Eberly technique by implementing the mapping 
li) |0) ^ li) (|0) + |3)) and other intermediate mappings 
of the form |i) |0) ^ (^1. li) + |T)) \^) §■ 




FIG. 1: Schematic energy level diagram for the combined 
harmonic oscillator/spin-1/2 system (only harmonic oscillator 
levels with n < 3 are shown). The arrows show the laser 
pulse sequence used to generate the state |J.) |0) from the state 
1*03) = U) (|0) + |3)) /^/2. The laser pul ses are applied in a 
stepwise fashion (labelled from 1 to 6) where the pulse areas 
i}ri,n't (marked next to the head of each arrow) are calculated 
according to the Rabi rate for the numbered transitions. In 
this notation, a "vr-pulse" would completely transfer all of 
the population from an initial state |i) in) to a final state 
It) \n'), for example. To generate I'I'os) from ||) |0), the pulse 
sequence shown here is applied in a time-reversed manner. 

The harmonic oscillator and auxiliary levels in our ex- 
periment correspond to the motional and internal states 
of a single ^Be"^ atomic ion trapped in a linear Paul 
trap ||l^ . We use the harmonic oscillator motional states 
along the trap axis (z direction) which are equally spaced 
in energy hy h x (2.9 MHz), where h is Planck's con- 
stant. In this direction, the ion is conhned by a static 
electric harmonic potential. The |t) and \l) (auxiliary) 
spin states are the F = 1, nip — — f and F = 2, mp = —2 



hyperfine levels of the ion's ^51/2 electronic ground 
state, which are separated in energy by approximately 
h X (1250 MHz). Applied laser radiation is used for state 
preparation and manipulation. A pair of laser beams 
detuned by approximately +80 GHz from the to 
electronic transition (A ~ 313 nm) drives coher- 
ent Raman transitions and couples the |t) and |J,) states 
and motional levels ||, Motion sensitive coupling 

is produced using non-coUinear beams with a wavevec- 
tor difference along z. The Raman laser beam fre- 
quency difference is tuned to drive It) \n — An) |t) |n) 
or It) \n) ^ It) \n) transitions, and the coherent transi- 
tion rate, or Rabi frequency, depends on both n and An 
|]l6| , |l9| . The experimental observable is the atomic spin 
state which we detect through state-dependent resonance 
fluorescence measurements at the end of every experi- 
ment §, |§. 

To demonstrate the Law/Eberly scheme, we 
configure the apparatus to generate the state 
1*03) = |i)(|0) + |3))/\/2 from |t) |0) using only 
transitions |t) |n — An) ^ |t) |n) where An alternates 
between and 1. The ion is initialized in the |t) |0) state 
with greater than 99.9% probability using stimulated 
Raman cooling and optical pumping pl} |. The six steps 
required to carry out the reverse process (produce |t) |0) 
from 1^*03)) are calculated according to the step- wise 
algorithm and are shown in Fig. 1. 

The state created after applying the Law/Eberly 
scheme is analyzed through measurements of Rabi os- 
cillations on the It) |n) |t) |n -|- An), An=0,±l transi- 
tions ||l|, The probability P| ^Y^il^nf to detect 
the ion in the |t) state is recorded after applying a laser 
pulse on one of these transitions for duration t. The ob- 
served oscillations (see Fig. 3) of P[ as a function of the 
laser pulse duration are fit to a sum of cosine functions 
with Rabi frequencies ^n,n+An{= ^n+An.71) constrained 
by the measured ratio of Rabi frequencies for the differ- 
ent motional levels ||2^ . The amplitude and phase (left as 
free parameters in the fit) of each frequency component 
are used to determine the probabilities |c|i | and |c|i | for 
i = 0, 1; 2, 3 iQ. We find that the observed ion popula- 
tion corresponds to the target state with 0.89 probability 
(Table 1), and that the populations in |t) |0) and |t) |3) 
are equal within the 0.03 measurement uncertainty. 

The probabilities |c|i|^ and |c|i|^ are also measured 
after each step in the procedure to generate 1^*03) and 
are compared to the theoretical predictions in Fig. 3. 
The Hilbert space trajectory from the initial to final 
state is somewhat complicated, with probability appear- 
ing, at least temporarily, in the |t) |n = 0, 1, 2, 3) and 
It) |n = 0, 1,2) states. 

The Rabi oscillation diagnostic determines the 
populations |c|i|^ and |c|i|^, but gives no informa- 
tion about the phase relation between the states. 
For example, measuring the populations in this 
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FIG. 2: Measured Rabi oscillations for the target state 
l^'os}. The probability to measure the atom in the \l) state 
is determined after applying a laser pulse coupling states 
II) |n) <-> It) |n + 1) for a variable length of time. Each data 
point (solid circles) represents the average of 600 experiments. 
The solid line is a fit to the data that is used to determine the 
populations in the first four motional states and the two spin 
states. Typically the fit determines that the 1/e time constant 
for the exponentially decaying envelope included in the fit cor- 
responds to 9 oscillations for the |i.) |0) ^ IT) |1) transition. 
The observed beating arises primarily from the oscillations of 
population in the \l} |0) and ||) |3) states, which have Rabi 
frequencies such that Qs^/Qqi = 0.60. The uncertainty in 
the spin state discrimination is smaller than the scatter in 
the data, which is mainly due to laser intensity and magnetic 
field fluctuations. 
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TABLE I: Measured state populations for the experiment 
with the target state l^'oa). Data similar to that in Fig. 2 
are used to determine the probability P = Icm^.^l^ to find 
the ion in the motional level n and the spin state nis =i,T 
for the intended target state l^&oa). This table shows the av- 
erage of populations determined from using Rabi oscillation 
measurements employing couplings with An = 0,±1. The 
uncertainty in the measured probabilities is 0.03 and is dom- 
inated by scatter in the Rabi oscillation data and the finite 
observation time. 



way cannot distinguish between the pure (coher- 
ent superposition) state described by the density 
matrix P = (U> |0) + |i) |3)) (( |( i | + ( 3 |( j |) /2 

and the mixed (incoherent) state described by 
P-(U>|0>(0|(i|+U)|3)(3|(i|)/2. To verify 

that our implementation of the Law/Eberly scheme 
establishes coherence we have performed a test ex- 
periment starting from \l) |0) using the target state 
I^-t) = 0.64 li) |0) -I- 0.77 It) |2), which would give 
kiol^ = 0.41 and Icpl'^ — 0.59. The first five pulses of 
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FIG. 3: Initial to final state Hilbert space trajectory. The real 
amplitudes ci„ and C|„ are shown for the sequence used to 
generate l^os). The solid lines are the theoretical prediction, 
while the data (solid circles) are derived from the probabilities 
|c|ip and |c|ip measured using the Rabi oscillation technique 
after each step. The amplitude was determined by taking 
the square root of the measured probability and assigning a 
sign consistent with the final state and the length of the laser 
pulses. 



the l^'os) sequence were used to generate \^t)- The 
measured probabilities for the experimentally generated 
state were |c|o|^ = 0.39, Icpl^ = 0.55, |c|o|^ = 0.03, and 
0.03 distributed among the remaining states. 

A coherent analysis pulse was applied on 
It) |ri -t- 2) ^ It) |n) transitions after the \'^t) state 
generation pulses but before spin state detection. The 
laser pulse area was adjusted to be a "7r/2" -pulse for 
the It) |2) ^ It) |0) transition. For the pure I^-t) state, 
the population would almost fully oscillate between the 
states It) |0) and |t) |2) as the phase of the analysis pulse 
was varied relative to the state generation pulses. No 
sensitivity to this phase would be observed if the state 
we generated was an incoherent mixture of populations 
(dashed line in Fig. 4). The measured probability to 
find the atom in the state |t) as the laser phase is swept 
is shown in Fig. 4. The amplitude of these oscillations 
can be related to the fidelity 

F = (*t|/5|*t) 

= 0.41/9^010 + 0.59p|2T2 + 0.495 {pioi2 + p-\2io) 

where p is the experimentally measured density matrix. 
We determine that F = 0.93 ± 0.03 using the measured 
populations and oscillation contrast to determine the rel- 



evant elements of the density matrix p as in ref 1 25 
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FIG. 4: Coherence fringes. After generation pulses for the 
target state 0.64 |i) |0) + 0.77 |t) |2) are implemented, a "7r/2" 
analysis laser pulse is applied with a controlled phase relative 
to the state generation pulses. Shown here are the resulting 
oscillations in the probability P| to find the atom in the \l) 
state as the laser pulse phase is swept. The amplitude of the 
oscillation determined from a fit to a cosine function (solid 
line) is used to establish the fidelity of the experimentally 
generated state. The oscillation centered around P| = 0.46 is 
consistent with the measured probability in the |t) |0) state 
(which is unaff'ected by the analysis laser coupling) and exper- 
imental error in the analysis laser pulse duration. The dashed 
line indicates the result if the prepared state is an incoherent 
mixture. 



In summary, we have demonstrated experimentally the 
scheme of Law and Eberly |^ for generation of arbitrary 
harmonic-oscillator states and its extension Q to arbi- 
trary harmonic-oscillator/spin states. The method can 
be generalized to higher dimensions [Q, to the genera- 
tion of arbitrary density matrices of harmonic oscillators 
[14 1, to the creation of arbitrary motional observables 
[15 1 such as the phase and to the generation of ar- 
bitrary Zeeman state superpositions P7| . The precision 
with which we can implement this technique has a direct 
relation to the efficiency of quantum-information process- 
ing using trapped ions . With sufhcient improvements 
in the fidelity of such operations, one can contemplate us- 
ing additional motional modes of motion as information 
carriers in this scheme. Of course, the same techniques 
can be applied in cavity-QED, the system in which it 
was originally conceived More generally, such tech- 
niques increase the variety of tools available for quantum- 
information processing and may eventually find applica- 
tion in areas not currently anticipated. 
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the manuscript. This work was supported by the NSA 
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tribution of the U.S. Government: not subject to U.S. 
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